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Phenomena (Kaiser-Wilhelm-Institut fur Stromungsforschung im 
Auftrage der Aerodynami ache Versuchsanstalt Gottingen, Institut 
fur instationSre Vorg&nge ) • 


I. INTRODUCTION 


The present treatise deals with the plane problem of the vibrating 
airfoil in supersonic flow within the scope of a linearized theory. 

The method is based on the application of results pertaining to the 
field of moving sound sources 1 2 and, in addition, on the introduction 
of the acceleration potential for the movement of the air particles. 

In view of the latter the method can be regarded as a direct extension 
of Prandtl's 'General Considerations on the Flow of Compressible 
Fluids '2 to the unsteady problem of the vibrating wing; the following 
makes use of all fundamental ideas of this theory. This extension 
leads to an integro -differential equation for the strength of the 
doublet distribution which represents the vibrating airfoil; this 
equation which is characteristic of the mathematical development of 
the theory, is solved by means of a Laplace transformation. The use 
of improper functions (first and second derivatives of discontinuous 
functions) which has become familiar in the modern quantum theory, may 
appear to be an innovation, but the fundamental idea which is expressed 
in this formal approach is nothing but the interpretation known from 
the theory of characteristics of hyperbolic differential equations, 
namely that disturbances and therefore also discontinuities in the 
hyperbolic case are propagated along real characteristics, in the 
present case notably on the boundary planes of the Mach wedge. 
Altogether the author believes to have found here a procedure which 
leads very simply and directly to a solution in the case of the 
vibrating airfoil moved in translation at supersonic speed; it can 
probably be profitably applied to similar problems as well. 

Dr. L. Schwarz has solved the present problem in "another way; I 
am particularly indebted to him for several critical observations 
which benefited my work. Apart from differences in the notation 
there is complete agreement between Schwarz 1 and my results so that 
their correctness ought to be fairly certain. I should like to mention 
that the theory has been carried through only to where it merges into 
the general final formulas of Schwarz. All further statements would 
be identical with Schwarz ’ statements and could not offer anything new. 


1 H. Hcnl , Ann.d.Phys. (now being printed) (1). 

2 L. Prandtl, Lufo 13 (1936) p. 313 I.2J TM 805. 
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II. STATEMENT OF THE PBOBIEM 

3 

As first noted "by Lanchester and carried further "by L. Prandtl, 

the introduction of an acceleration potential from which the 

— * ) 

acceleration vector b 


+ Ck, vsO 


(1) 


« 


Ql is used for German script 


of -tie air particles (h velocity vector) Is derived according 
to d * grad 1 }/ j proves to he especially useful for the treatment 
of the wing in a compressible medium. \{/ then results according to 
Euler 'b hydrodynamic equations for frictionless liquids, apart fr can. 
an arbitrary constant, as identical with the negative pressure 
function, thus 



a Constant 


( 2 ) 


The fundamental Idea in applying the acceleration potential to the 
wing theory in the shape given to it by Brand tl is the following: 

If the perturbation velocities are linearized, \|/, on the one hand, 
satisfies the well known equation of wave propagation 


A\}/ 


1 a 2 * 

7 57 


(3) 


(c sonic velocity) on the other hand, however, experiences 
according to ( 1 ) a discontinuous change in penetrating the lifting 
surface on which there appears a pressure jump between the lower and 
upper surfaces. These conditions suggest the use, in determining * 4 / , 
of analogies from electrostatics (for the steady-flow problem) or 
electrodynamics (for the unsteady-flow problem) . 

Wo will limit ourselves to the two -dimensional problem of the 
infinitely long wing and suppose that the lifting surface is a strip 
of width l, located in the x,z -plane, parallel to the z-axis. 

Furthermore, it will be useful to employ a system of coordinates 
which- is attached to the airfoil so that the wing in reference to 


■S 

L. Prandtl, elsewhere. 
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this coordinate system executes small vibrations but has no trans- 
lation; if one assumes that the wing in the original coordinate 
system (medium at rest) moves along the negative x-axis with 
the velocity U, the wing in the new coordinate system (without 
translation) is in a flow of the same velocity in the positive x-axia, 
and equation ( 3 ) must be replaced by 

-° " (3,) 


In this equation, d/5t is to be left out in the case of steady 
flow. For the problems of unsteady flow we will limit ourselves 
to the case of the harmonically vibrating (infinitely thin) wing 
and accordingly make the substitution 


^ * <p(x,y) a -la)t 

The wave equation ( 3 ') changes into 


<*) 


(*f - 1 ) 


^ 9 - + £ifcM 




I) 9 - 0 


( 5 ) 


with 


k = 


o> 


m-2. 

c 


(5a) 


(k quantity propagated^ Mach number) 

We now assume that $(x,y) is a solution of the wave equation ( 5 ) 
we call it an elementary solution - which corresponds to a sonic p«int 
source located at x * y = 0 .* 


4 


H. Eonl, elsewhere. 



1 


NACA TM 1238 


5 




Then obviously the "source potential" 

« U f l $(x - |,y) f(£) dg 

Jo 


A * 

with an arbitrary source distribution U f(x) along the wing 
chord 0 £ x £ l anywhere outside of the lifting surface will 
satisfy the equation (5) and the seme is valid for 


,.*9b 

5y 


& $y(* - §»y) f(g> *£ 


.0 


( 6 ) 


(The faotor U is added to avoid unnecessary factors in the 

* following formulas) . The latter expression obviously represents 

apart from the time factor e'i^t -the retarded potential of a 
doublet distribution; f (c) ie now a measure for the strength of the 
doublet distribution and is, therefore, proportional to the Jump 
in <p while penetrating the surface. If we substitute (6) into ( 4 ), 
will meet all requirements which are to be demanded of the 
acceleration potential and can therefore be identified with it. 
However, it muBt be stressed that \| / and <p will be uniquely 
determined for the whole space only when the source and vortex 
distribution is given everywhere and in addition a boundary condition 
or ray condition ® at infinity is prescribed for <p. In fact, a 
vortex wake in general adjoins the trailing edge of the wing which 
also will contribute to the acceleration potential. The. conditions 
for the case of supersonic flow are rendered especially simple by 
the fact that each source or vortex element takes effect 
only and that therefore the air forces on the wing are solely 
determined by the doublet distribution replacing the pressure jump. 
Concerning the boundary condition at infinity we shall have to require 
that <p and therefore the perturbation velocity components u and v 
disappear outside of idle Mach wedge which starts at the leading edge 
of the wing (therefore, especially for x<0) but in the in finite 
inner region of the Mach wedge satisfy a ray condition; these 
requirements are automatically met by our equation ( 6 ) according to 
the definition of the elementary solution. 

t Therefore, our task will be to find conditions which uniquely 

determine the source distribution t(£) in (6). One such condition 

S' " ‘ — ‘ '■ ■ ■■ ■ — — — ... 

Compare A. Somnerfeld, Jahreeber. d. EMV 21 , pp 309-353 
0 - 913 ) [ 4 ] • 


\ 
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results from the fact that the dovnvash distribution 

w(x, t) * v 0 (x) e itUffc on the wing is given. Therefore, ve further 

need at first relations between the acceleration potential t and 
•the perturbation velocities u, v (u horizontal, v perpendicular) 
which permit the expression of the latter by within the scope 
of the linearized theory 6 . Thus we obtain corresponding to our 
assumption about the establishment of the coordinate system according 
to (1) apart from terms of the second order in the perturbation 
velocities 

(7) 

3t Sx 


and, therefore, by substituting 1? = grad t , v = grad ^ 

sv a® 

+0C(t) (7') 

0*b OX 


where ^ (t) stands for an arbitrary function of time. Since we 
will assume, however, that no change shall take place at infinity, 
we can immediately specialize 'X(t) » C. If we assume t (x,yj t) 
as given there follows now for the velocity potential 7 

t - 3£ r iL ') ax ' + f i + # .. {8) 

with D as a further (insignificant) integration constant. In this 
equation the integral appearing in the first term 



6 Compare H.G. KGssner, Lufo 17, p. 370, Abschn.4 

7 An arbitrary function f (t x/U) to be added in (8) according 

to (7') is reduced to a constant since otherwise the wave equation 
for *W (analogous to (3) and (5), respectively) would not be 
satisfied. 
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corresponds to the potential of the perturbation velocities, the 
^ part containing x linearily corresponds to a superposed undisturbed 

flow in the x -direction; therefore one has to equate C * U*. 

From (8*) follows 


8 ? 
V » — 

by 



(9) 


as well as a further integral for u not needed at present. 

If we make use of (4) for harmonic vibrations we obtain for 
the x-axie: 

/ V -lot a "lot 

v = v_(x) e = Limes 2 

O' TT 

y — » o y 


or 

e ^ v Q (x) = Limes -i- f 9y(x’,y) e dx’ 

y — *o v <r<» 


[ 9 y W,v) e^^^dx 1 (10) 

v '*cc 


and finally by differentiation with respect to x: 


e 




= Limes ~ <P (x,y) 
y — 40 u J 


( 11 ) 


We now have to substitute (6) on the right side of (11) . Taking the 
fact into consideration that only those s'ources contribute to the 
downwaah at the point x of the wing in supersonic flow which lie 
upstream with respect to x, for which therefore - x, we can 
replace the fixed upper limit of the integral (6) by the variable x; 
as the lower limit we will select -w, establishing f(£) = 0 
for < 0, so that f(£) for ^ = 0 will, in general, be discontinuous. 
In the limit y> — >0 the right side of (il) must approach the same 
limit value for both signs of t) * jyj. Thus we obtain from (6) 
and (11), taking (5a) also into consideration, 

I 0 w (x ii) f(£> <£ - S 2 - 1 fi T 0 0*) . 
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The right side of (12) is a given function of x; however, 
it must, be taken into consideration that now v Q (x) also at the 
point x * 0 is to be regarded as discontinuous (v 0 (x) = 0 

for x<0), so that at x a 0 has to be represented by an 

dx - - 

improper function (5 -function) . Correspondingly, in the deriva- 
tive £>yy also singularities of higher order appear in the integrand 

By the integral equation (12) with a given kernel $ yy and given 

right side the distribution function f(x) is uniquely determined- 

The calculation of the pressure Jump on the wing which is the 
real goal of this paper, results after solution of the equation (12) 
immediately from (2). First, the change of the pressure function 
corresponding to the infinitesimal pressure change Ap will be 



with p 0 standing for the average density of the medium, and 

therefore, according to (2) the pressure Jump n between lower and 
upper side of the wing (the indices + and - refer to upper and lower 
side) : 
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* * Ap_ ” Ap + = P 0 (vj/ + - t.) 



£ 

• Equation (13 ) for the pressure Jump « is essentially 
identical with Bernoulli 's relation between pressure and velocity j 
this relation is represented within the scope of the linearized 
theory for the case of unsteady flow by the equation 


/ 


P Sx 


at 


= Constant 


(or) 


or, using the notation of 


u 


M. 

5x ’ 


by 





= U(u + - u.) 



CP) 


The calculation of the horizontal component u yields according to (8) 
at first 


u 


b'ir 

5x 


= -^-'Kx,y;t) 





dx’ 




x,y;t) 


8t 


(r ) 


therefore, according to ($) in agreement with (13) 


u 0 **+-+- 

The explicit calculation of u and . from f according to (3) 

8t 

proves therefore to be a detour here. 
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Thus the solution of our problem is traced back to the solution 
of the integral equation (12) for the’ distribution function f(x) 
from which it can be calculated according to (4), (6), and (13)* 
Since singularities which originate in discontinuities of v*(x) 
and $(x, y) appear on the right as well as on the left side 
of (12), the solution of (12) shall be preceded in the following 
section by a discussion of the properties of the 5 -function. 


HI. THE METHOD OF THE 6 -FUNCTION 


By the function 6(x) - also written 5(x^y) - ve- understand 
an improper function, which, with the exception of the points x =* 0 
becomes infinite in such a way that 

AXg 

1 5(x) dx « 1 (14) 

• 1 

becomes valid if the interval . . . Xg contains the point x * 0 

as an inner point. In order to attach a stricter mathematical meaning 
to the arithmetical operations making use of 6, it is necessary to 
consider 6(x) as limit of a sequence of continuous functions. We 
assume -therefore at first continuous positive functions ?p(x) which 

have a steep maximum at the point x = 0 and are also subjected 
to the condition 

<p p (x) dx * 1 (14 ' ) 


Let <Pp(x) be a function of a parameter p in such a way that <p p (x) 

everywhere, with the exception of the point x » 0, approaches zero 
as p approaches a fixed limiting value, for instance p «* 0. The 
"limit function" 5(x) thus characterized may be denoted as "prong 
function." All mathematical expressions containing the 5 -functions 
in a symbolic way must, therefore, be considered as limiting processes 
in which at first 5(x) is replaced by the sequence of functions <p p (x) 

and finally the limit p— *o is taken. 



i/x 
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After this c cement there results immediately the property of 
the 8 -function which is most important for its application, the 
integral relation ■ 

f*2 

I 5(x - a) p(x) dx * F(a) (15) 

for x^ < a < Xg, in case F(x) at the point x = a is continuous. 

In addition, we will consider the sequence of the deriva- 
tives 9^(x) of the function $>^(x) defined above. We call the 

limit function characterized by this sequence 8^(x); therefore, 
we may equate symbolically: 




(*> 


d8(x) 

dx 


(16) 


Since we may transform integrals, which contain in the integrand 
the continuous functions by integration by parts into 

integrals which contain q>p(x) instead of <p ’(x), there follows 

from (15) also 



- a) F(x) dx = F’(a) 


(IT) 


■ where < a < Xg, for all functions F(x) which have a continuous 

derivative at the point x * a; use is made of the fact that the 
functions <p (x) disappear in the limit p— * 0 at the limits of 
the interval. 


Finally we mention that the sequence of functions <p^(x) 

also be regarded as a sequence of the derivatives of function 
which approximates the "step function" 


may 

gp(x) 


0 for x< 0 
for x £ 0 


s(x) - <j 

I 


1 


(18) 
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In a continuous way. In this sanse we may then substitute symbolically, 
in analogy to (11 }, 


»(.) = - (») 


The application of the method of the 5 -function to the present 
problem is based on the fact that the elementary solution <5 (x, y) 
of the wave equation (5) is, for supersonic conditions, quite generally 
a discontinuous function - namely zero outside of the Mach wedge, 
but fO inside of it and also on the boundary when approached from 
the inside so that the first and second derivatives <$ y (x, y) 

and 0 yy (x, y) which appear under the integral according to (6) 

and (12) lead to the separation of improper function parte on the 
boundaries of the Mach wedge. On the other hand, the performance 
of the prescribed integrations according to ( 15 ) and (1?) does not 
present any difficulties as long as we may suppose the distribution 
function f(x) and its derivative f *(x) everywhere as continuous . 

This assumption' is fulfilled for the extent of the wing 0 < x < l 
with exception of its leading and trailing edge as well as in front 
of and behind the wing, that is, for x < 0 and x > 2, where f(x) * 0. 
The points x * 0 and x * l are exceptions. We will try to evade 
this difficulty by assuming at first that f(x) is continuous also 
near the points x = 0 and x = l and will proceed only subsequently 
to the limiting case of a discontinuous distribution function. This 
procedure always leads to a uniquely determined result. The mathe- 
matical meaning is that one has to proceed first to the limit p — » 0 
and only subsequently to the limiting case of a discontinuous 
distribution function or downwash distribution on the wing (right side 
of (12)) since the performance of the Integrations (6) and (12) already 
involves the application of the 8 -symbol. Such a procedure that regards 
the compression shock along the boundary of the Mach wedge which 
corresponds to the elementary solution $’(x, y) as infinitely sharp 
(linearized theoryl) in contrast with the always incompletely realized 
sharp limitation of the wing, is probably fully Justified from the 
physical point of view. 


17. THE PLATE AT A GIVEN ANGLE OF ATTACK 
IN SUPERSONIC FLOW (STEADY~FL0W PROBLEM) 


We illustrate the method for calculation of the air forces, at 
first by the steady flow - example of the plate at rest at an 
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infinites Inal angle 3. We shall obtain again results already 
known. 9 

In this example the dovnwash is constant at "both sides of 
the plate, (we equate v Q = - b), but zero in front of and behind 

the plate. Obviously 3 =— (apart from quantities of a higher 

u dv 

order) (See fig. 1.) If we now form the derivative , as provided 

oJC 

for in (12), this quantity will become at x = 0 and x « l, 
respectively, -«® and + », but in such a way that the function 

f x av ! 0 for x < 0 

v Q (x) . J Lois, dg, » < -b for 0 < x < i (20) 

d £ L 0 for X >l 


at x » 0 and x = l experiences a discontinuity by +b. According 
to (19) we, therefore, may substitute for the right side of (12) - 
taking into consideration that the second term of the right Bide 
disappears in the steady-flow case (oj = 0 or k = 0, respectively) - 


d v Q 




ss 


-b5(x - 0) + b8(x - l) 


(20') 


On the other hand, let us consider the elementary solution $ (x, y) 
appearing in (12) and its derivatives. By <±>(x, y) we understand 
in the steady-flow case a discontinuous function which assumes inside 
of the Mach wedge 


y 

x 


tan a, 


1 

x > 0 : tan a = ' 7-A- — ' 

Vvr - 1 


( 21 ) 


with the apex x = y = 0 and at half the opening angle a a constant 
value (that is, a value independent of x, y and with a certain 
normalization dependent only on M) , but disappears outside the 
Mach wedge. 10 It does not mean a restriction of universal validity 

Compare Ackeret, Zs.f. Flugtechn. Motorluf tsch . Bd. 16 (1925), 

P.72 [61 

"’Compare L. Prandtl, elsewhere. Equation (30). 


Ik 


NACA TM 1238 


for our problem if we equate the value of the function inside the 
Mach wedge - 1. Introducing the step function (l8) obviously 
makes 


f s(x + cot a y) for x > 0 
$(x,y) » < 

jo For x < 0 


(28) 


with the upper sign valid for y > 0, the lower sign valid for y < 0. 
If we further take (19) and (l6) into consideration, we obtain 
for x >0 with the abbreviation q ® x jp 00 1 a y: 


$y(x,y) « + cot a $'(<l) ® + cot a 8(q) 


(23a) 


0 (x,y) c cot^aO"(q) * cot^a B^^(q) (23b) 

yy 


If we now apply the integral relation (17) the integral on the left 
side of (12) becomes, if we use again n * |T| instead of y. 


Px 


©yy(x - g, *tj) f (g) d£ - cot 2 a 


d-oj 


nx 


lj-0B 


5 (1) (x - g - cot a q) f(g) d£ 


* oova f '(x - cot a q) 


( 2 k) 


Therewith we obtain from (12) with the help of (20’) in the 
limit t) — 3.0 the differential equation 

cot 2 a f '(x) = - b8(x " 0) + b6(x - l) 


(25) 
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from which we conclude that f (x) must he discontinuous at the 
points x * 0 and x * l; and so, apart from a constant, the 
function 


f(x) * as(x - 0) - as(x - l) 


(26) 


is, according to (19), the solution of (25), if a is selected: 


a 


-h 

cot a 


.(26*) 


The constant, which is as yet undetermined, is obtained from the >• 
fact that f(x) must always vanish ahead of the leading edge (x = 0) 
and is, therefore, left out in (26). Of course) we likewise 
put f(x) * 0 behind the trailing edge (x > l) , a condition which 
is then automatically fulfilled according to (27) . However, there 
is obviously only one integration constant at our disposition for 
the transition from (25) to (26). The reason for the possibility 
of giving the solution of (25) in the form (26) is, according to (20), 
that the dovnvash v Q performs an equal and opposite jump at x = 0 
and x = l so that, behind the trailing edge v c = 0 . Therewith, 
we only confirm the correctness of our original assumption- For 
non-plane lifting surfaces or for 'the problem of unsteady flow the 
above statements will, in general, be no longer valid; a vortex **»ke 
will then adjoin the trailing edge 

We calculate the horizontal component u of the perturbation 
velocity- From (8*) and (4) one obtains for the steady- flow case 


u 


Sx ' V 


ir v(*,y) 


and further, according to (6), 


u 


ix 

V x "■£» d £ 


( 2 ?) 


(27a) 


u 
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if ve fix the integration limits as in (12). If we go back to (22) 
and (23a) and use equation (15) for the transformation of the integral, 
we obtain, because of (26) and (26a), 


u * *cot a 8(x - g - cot aq) f(§) 


* $cot af(x - cot an) 

a ± — ^ — s(x - cot an) * s(x - l - <?ot an) 

cot a L 


and in the limit tj — » 0 : 


u 0 « t — 7 — 
0 cot a 


s(x * 0) 



If we compare (28) with (20) for which we can also write 

« -b ^e(x - 0) • s(x - l) 

we find 


(28) 


(20a) 


Uq • |V Q tan a 


(29) 


(In all formulae the upper sign always refers to the upper side, thp 
lower one to the lower side of the lifting surface . ) 

The pressure difference n on both sides of the lifting surface 
is now immediately obtained from (13), (27), and (29). There follows 
with v Q «= -b: 


rt * 2p Q Ub tan a 


(30) 
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and from this equation the lift 


A s nl s 2p 0 UbI tan a 


(30a) 


In order to obtain the lift coefficient c a we hare to divide 

by -£pU^l. If we replace b by the angle of attack 0 « , 
there follows from (30a) 

c = 40 tan a (31a) 

The result for the drag coefficient c„ will be, since the auction 
force at the leading edge is eliminated In supersonic flow, 

c v = 0c a - 40^ tan a (31b) 


The formulae (31a) and (31b) are in agreement with the results of 
Ackeret who for the first time calculated the air forces on the plate 
at a given angle of attack in supersonic flow in 1925. 

As a conclusion to this treatment we want to stress once more 
that the velocity and pressuro distribution of the plate at a given 
angle of attack were traced back essentially to the die contlnui ties 
of the distribution function f (x) . This fact is expressed clearly 
in the connection between the formulas (26)and (28) . Thus a discon- 
tinuity surface of the u-distribution as well as the v -distribution 
spreads from the leading edge of the lifting surface which coincides 
with the boundary planes of the Mach wedge and causes an expansion* 
wave in the upper region y > 0 and a compression wave in the lower 
region y <0 (fig- 2). The effect of this discontinuous change 
in velocity is in our example exactly cancelled by the second discon- 
tinuity surface starting from the trailing edge x = l (compression 
wave for y > 0, expansion wave for y < 0) so that the flow behind 
the second Mach wedge which starts at the trailing edge is again 
undisturbed. Since obviously the sphere of action of a discontinuity 
surface of f (x) in supersonic flow extends downstream only, there 
follows that the problem of pressure distribution on the wing in 
general is already uniquely determined without ea v additional "flow-off 
condition" (requirements for finitenese for the v-distribution 
immediately behind the trailing edge). The discontinuity of f(x) 
at x » 1 involves in general the appearance of vortices at the 
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trailing edge which, however, need not he ascertained for the 
investigation of the air forces on the wing. 


V. BENDING OSCILLATIONS OF A PLANE WING AT ZERO ANGUS OF ATTACK 
SETTING UP OF AN INTEGRO -DIFFERENTIAL EQUATION 


The tending oscillations of a plane wing at zero angle of attack 
represent the simplest problem of unsteady flow. Let y Q be the 
(infinitesimal) displacement of the mean camber line from its zero 
position (x-Axie); we then substitute 


v m Ae for 0 = x = l 

0 


(32) 


Then the downwash in the same interval becomes 


v(x, t) * v Q (x) a = -icuAe ~ iayfc 


and, therefore, 


^1° a -i<jjA5(x - 0) + • . ■ (33) 

By ... we refer to the contribution of the second discontinuity 
of v at x = l which, however, does not explicitly enter our 
calculations since we can break off the wing at an arbitrary point 
x = l without having to change the results for x< l. 

Our most important problem now is the transformation of the 
integral on the left side of equation (12) and the transformation 
of this equation into an integro -differential equation for f of 
the ordinary type. The integral equation (12) has as in the problem 
of steady flow an essentially degenerate kernel since the elementary 
solution $ on the limit of the Mach wedge becomes discontinuous by 
itself; its derivatives, therefore, are to be represented by the 
improper functions 8 and 8 (1), Therefor© we separate 0 (x, y) 

and $ (x, y) each into a degenerate part due to the discontinuity 

v if 

of the function 0 at the boundary of the Mach wedge and into a 
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regular part due to. the analytic behavior of 5 in the inner region 
of the Mach wedge. The elementary solution $ (x, y) of (5) for the 
time -periodic case of a sonic point source can he taken from a 
former investigation. 11 We normalize the function $ in such a way 
that the Jump of $ in penetrating the Mach wedge has exactly tho 
amount 1; the elementary solution then reads: 


C (*,y) 





inside the 
Mach wedge 


0 


outside the 
Mach wedge 


(3*0 


where we further introduce for abbreviation 


€ 


2 

cot a 




1 


(35) 


Taking the discontinuity of ft on the boundary planen x - cot ajyj = 0 
of the Mach wedge (for x > 0) into cons ide ration, we obtain 


$ y (x,y) = +,/s 5 (x T\j e 1 ^ + fc*(x,y) (36a) 

® ( x >y) - (x -\/€ n) +0* (x,y) (36b) 

yy 


where § *(x, y) characterizes the part of $ (x,y) which is, within 
the Mach wedge, continuous (and analytic); outside of the Mach wedge 
the derivatives disappear. With the insertion of (36b) into the 
integral on the left side of (12), the first part of (36b) yields 
by integration by parts 


11 

H. H8nl, elsewhere. Equation (4o)' 
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OX 


1-03 


8 (l) (x - £ - iff i|) e&fr-g) 


*(£) <*£ 


nx 


+c 


S(x - £ - & 1 \) 


d SL 




f(£) 


0-« _* 

{[- i & f(g) + f *(£)] 


£=x-v?r| 


and therefore in the Halt n — »0: 


- lMkf(x) + € f '(x) 


The integral equation (12) therefore assumes, on account of (33), 
the form 


(37) 


The right aide of (37) shove that f(x) in any case must he dis- 
continuous at x 0 0, ao that f '(x) contains an Improper function 
8(x - 0) as a constituent part that with proper choice of a 
multiplying constant will cancel out against the right side. This 
is natural since the function f(x) for x < 0 , that is, ahead of 

the wing, must vanish together with its derivative f *(x) (so that 
(37) for x <0 is solved hy f(x) » 0), and therefore will, in 
general, show unsteady behavior at the front boundary of the wing, 
as already shown in the example of steady flow (IV). If we, therefore, 
set up the equation 


- iMkf(x) + ’(x) + 1 ^ (x - £ ,0) f(£) <3£ 


J-00 


j'- ioA 6(x - 0) - i ~^j for x » 0 


for x <0 


f(x) * as(x - 0) + f^(x) 


( 38 ) 


where a represents a constant yet to be determined and the 
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function f x (x) defined for x ^ 0 has to satisfy the condition 
that it changes continuously into 


f x (0) = 0 


( 3 &) 


(37) 150001008: 


- iMka - iMkf 1 (x) + eje&(x - 0) + f^(x)] 

+ a f X 5yy( X * *£ + r ®yyCx - g ,0) f x (£) d g 

Jo Jo 


-itflA 


S(x - 0) - i ~ 


( 39 ) 


If we now put 


a 



(to) 


the 8 -function will he completely eliminated from (39) and we 
shall obtain for fi(x) a linear Inhomogeneous (Volterra) 

Integra -different lal oqy^frlga <?,£ convolution type, pf.ga 

otherwise, regular character 


/ P x 

lMkf^^Cx) + Cf^x) + $*y(x - g ,0) f x (^) d£ 

Jo 


I* 


a( iMk / ®* y (£,0) ig I “ if ca 




(to) 


its solution is uniquely determined by the boundary condition (36a)* 

The solution of the Integra -differential equation (to) 
required a further determination of its kernel. If we go back to 
(3C) we obtain by differentiation for the inner region of the 
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Mach wedge: 


0*(x,0) * 0 


^yy( x >0) 


€ 


J lfe) 

kx 


(42) 


where the equation 


Jq(z) = *J x (z) 


was made use of. If we introduce further the dimensionless variable 
(reduced length) 


z 



(43) 


instead of x, we can also write: 


fc* (x,0) G(z) ; G(z) * iilL (42a) 

z 


Equation (4l) then assumes the form 


-iMF^z) + P^(z) + 


Pz 


\p 


0(z -£) F x (g) d £ 


a 


iM 



G<0 d t\ 


(44) 


The symbol F 1# rather than fj^was selected in this equation in view 
of the replacement by z of x. In the following section VI we are 
going to solve the lntegro -differential equation (44) for the 
boundary condition (38a) by means of a Laplace transformation; 
the result reads (for z * 0): 


v*> - * i '>“ z v z > r w - 1 


1 * 


(45) 



NACA TM 1238 


23 


With this- the distribution function f(x) also is determined. By 
-writing f(z) = aH(z) instead of f(z), we obtain from (38) 


F(z) 5 aE(z) * 


- 4 f 

Jo 


e iMu J ft (u) du 


for z ^ 0 


(46) 


for z < 0 


where the value (4o) has to be substituted for a. 


In addition we shall prove that the pressure distribution on 
the wing for the general case of unsteady flow also is determined by 
the distribution function f(x) or F(z). respectively, alone ■ 

If, to this end, we go bach to (4) and (6), we have at first 


\{/(x,yjt) =. e“ ioit U 


nx 


$ v ( x -S, ±r\) f(|) 


'•CO 


and from this equation according to (36a), because ox‘ the disappearance 
of $* y (x,y) in the limit jy| — *0, (compare (42)) 


ty(x,±0;t) = e” iayl ' Limes 

r\ — >0 


\x 

U-co 



jfVe f(x) e 


•last 


(47) 


By substituting (47) into (13) we thus obtain, again replacing 
f(x) by F(z), for the pressure jump it on the wing 

it ■ -2p 0 UVT F(z) e“ icBt (48) 

This formula is valid in general. From it we obtain for the special 
case of bending oscillations, according to (46) and (4o), 
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"by which the problem of the calculation of the pressure distribution 
on the wing is, for this special case, solved. 

In discussing the equations (46), (47), and (49) one must not 
overlook the fact that the variable z contains the quantity <, 
which decreases with M — >1 to zero, according to the equation 

W 

z » £ . Therefore, the graphs for the F-, u-, and n -distribution 
as functions of the chord x in the limiting case M — > 1 shrink more 
and more toward x • 0. 


VI. SOLUTION OF THE INTEGRO -DIFFERENTIAL EQUATION IN 
THE CASE OF AN ARBITRARY DOWNWASH DISTRIBUTION 
BY MEANS OF A LAPLACE TRANSFORMATION 


Before we turn to the general problem named in the title we 
shall first derive the solution for the special case of the bending 
oscillations by means of a Laplace transformation. The formalism 
of the calculation can then be transferred without any difficulty of 
the general case of an arbitrary time -periodic downvaeh distribution. 

^ Bending Oscillations as a Special Case 

Our problem is the solution of the integro-differential 
equation (44) for the boundary condition F, (0) =0 by means of 
a Laplace tronsf ormation . 

If we denote the Laplace transformation L g (F-]_(z£] by 
<p(s), L [G(z)J by g(s), the equation which corresponds to (44) 
is linear in <p(s) and becomes, taking into consideration the 
boundary condition F^o) = 0 13 


[s - iM + g(s)j <p(s) = aig 

iM - lg- - (“ O(J) 


- Jo - 



4 ~ 


(50) 


id 

It is noteworthy that L g [F'(z)J * s<p(e) » F-^(O) 

compare G* Doetech, Laplace transformation (J. Springer 1937), 
P- 153. 
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13 


Tho Laplace transform of J- L ( z, )/ t i s: 


^ («) 


L 2 J 


ICO 

iiM .' ez i* . \/?n 


and therefore, (compare (42a)) 


g(e) 


■JOO 

ilM e iMz e’ sz ds » ^(b - iM) 2 + 1 - (s - IM) (51) 


u° 


The solution of (50) for cp(s ) is, therefore. 


<p(s) <= a 


iM - i~ - g(s) b - \J(s - 1M) + 1 - i^ 


&s - iM + g(s)| s t/(s - iM)^ + 1 


(52) 


As is veil known, the inverse transformation cp(s) — » F 1 (z) is given 

hy 


F 1 (z) * l ” 1 [*(.)] 


1 

2iti 


e sz (p(s) ds 


(53) 


where the closed path of integration traversed in the positive sense 
encloses, in a unique sheet Tschlichten Blatt) of the complex 
s -plane, the singularities of the integrand located in the finite 
part of the plane. If we finally introduce t * s - IM as a new 
integration variable, F^(z) becomes, according to (52) and (53) 


F x (z) 


a e 
2«i 


1 e zt 


/ 


/?Ti t + m 


- <£ 


Hi 


(t + iM 


A 


dt 


14 


G. Doetsch, p. 313 and 403. 


( 54 ) 
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0ns obtains the integral of the first term of the sum immediately 

from the integral representation for the Bessel function J : 14 

o 




.tz 


2iti 


/FT 


dt 


(55) 


the integral of the second term of the sum can be evaluated 
immediately, according to the residue theorem; for the integral of 
the third term of the sum one obtains Imm ediately because of 


"Si' 


e lMz 

2ni 


/ ( t ' + tTTFTi 4t ] ‘ "S?T 


dt » 

(55a) 


the equation 
lMz 


»tz 


2«i J (t + iM) jt d + 1 


dt a I © :3Jviu J (u) du + C 


(55b) 


where c » 0 since the integral for z « 0 must vanish according 
to the residue theorem. Hence (54) becomes (45), as had to be 
demonstrated* 


{*>) General Case of the Downvash Distribution on the 
Harmonically Vibrating Wing 

In this section we assume an arbitrary down wash distribution 
of the form 


v(x,t) =-£- + v (x) e~ ia)t 

St 3x 


( 56 ) 


14 Courant and Hilbert, Methoden der mathematischen Physik I, 
(verlagj. Springer), p. 390. 
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with, hovever, the distinct requirement that the (infinitesimal) 
displacement y Q of the mean camber line from its zero position 

(x-axis) must everywhere he continuous and especially for x » 0 
must change continuously into y 0 «= 0. For x < 0, we assume 
y Q = 0, and therefore v 0 (x) = 0. It is easy to allow for 

discontinuities of hy passing to the limit of a sequence of 

continuous functions y Q (compare case c). 

For the general case of an arbitrary dovnwash distribution the 
right side of the integral equation need not be further specified; 
also, in the equations corresponding to ( 37 ) and ( 39 ). the same 
general expression appears on the right side as in (12). If we 

loc 

again introduce the dimensionless variable z » -== instead of x 
and denote v Q (x) as a function of z by y(z),ve obtain 


dv o(*) . ,k k Ay(z) 
dx “ e d z 




The form of the integro -differential equation for F(z) E f (x) 
which corresponds to (44) becomes, therefore, for positive values 
of z (after again dividing by a factor k) 

-iMF(z) + F *(z) + f G(z -£) F(£) d£ «-|r(z) (57) 

J° 

with "Hie abbreviation 

r ( z ) . Mil - i^y(z) (57a) 

Equation (57) has to b© solved subject to the boundary condition 
F(0) *» 0, as can be inferred from the assumed continuity of y_ 
at x * 0; for z < 0, F(z) =0. 

If we apply the Laplace transformation to ( 57 ) and denote 
the Laplace transform of F(z) by tJ/(s), we obtain 

[s - iM + g(sfl \j/(s) = [r(z)J 
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If we carry 
taking (51) 


out the inverse transformation \{/(s) — *F(z), we obtain, 
into consideration, according to definition 


1 


F(z) 



- 1 < 


1 

|/(s - iM) S + 1 



( 58 ) 


If we now consider that, according to ( 53 ) and ( 55 ) 



- 1 ■ 

l/(s - iM) S + 1 


= e^ J o (z)£ 5 X(z) 


( 59 ) 


and furthermore, that the result of the inverse transformation of the 
product of two functions (in s) is the convolution of the inversely- 
transformed factors (In z), thore immediately follows from (58) 
and (59) 


F(z) 



iM(z-£) 


J (z 
o 


t ) r(J) d£ 


(60) 


which represents the solution of equation ( 57 )* By inserting ( 60 ) 
into ( 48 ) one then obtains for the pressure Jump on the wing 
(for z = 0) 


it 


•SPo 11 


_1 

ft 



K(z - £) r(£) 



(61) 


with the meaning of K(z) explained by (59) 


(c^ Transformation of the Results (60) and (6l) 

We at first confirm that the special formulas ( 46 ) and ( 49 ) for 
bending oscillations as a limiting case are contained in our general . 
results ( 60 ) and (6l). We have to make allowance for y(z) Jumping 
near z = 0 practically discontinuously from 0 to the constant 

dy(C ) 

final value a€ = -iaA so that in (60) the term : in r(£), 

equation (57a), makes a contribution only at the point £= 0 and there 
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can tie represented by 


K 


a* 8( i - 0) 


Thus we obtain according to (15) 


K(s - £) MSI » ac I K(z - £) S(£ - 0) d£ *= a € K(z) 


F(z) therefore becomes, according to (57 a ) and (6o), in agreement 
with (46) 


n ! 

1 K(z - £) d£ = a 

, pz 

K(z) - K(u) d 

Jo 

^ fcO J 


In order to recognize that the general case also can be traced 
back to the type of bending oscillations we perform an integration 
by parts on (60 ), introducing temporarily the function 


in 


K(z) = / e iMu J (u) du (for z = 0) 


(62) 


Thus we find 
*z 


K ( 2 * £) 7(g) <*£ *'lC(z *£)>($) 


♦ 1 Z k(z - t ) d l 




The part to be evaluated between the limits vanished, since 7(0) = 0 
and K(0) = 0. Therefore, (60) may be transformed according to (57a) 


r(*) -7 


K (I - t ) -4 k(> - £ )j-^d£ 
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and from (46) and (62): 


P(») 





(63) 


Only the influence function H(z), which is characteristic for the 
"bending oscillations, appears in (63) besides dy/d£. Thus the 
pressure distribution appears composed of the integrated effect 
of infinitesimal steps, the magnitude of which is proportional to 
the derivative dy/dz . 

For practical needs it will, in general, be more useful to 
eliminate the derivative dy/d£ in (60) by integration by parte 
and to trace the pressure distribution back to an integral with respect 
to r(0* We find 

lz 

K'H) HO K 


since 7(0) » 0 and K(0) » 1. If we take the derivative of K we 
obtain from (59) and (60) after an easy transformation for z? 0 

I) 

»(«> ■= T ’'<*> + 7 I [i * 0 (‘ -£) -Ji( s -t>] «£> d £ 

Jo 


K(z 




y(z) + 


In this formula there appears nothing but the downwash distribution 
y(z) and a characteristic influence function. • 

The final formulas (60) to (64) are in complete agreement with 
the results of Mr. L. Schwarz 15 , as one can easily confirm after an 
adequate change of the symbols. Therefore, further detailed statements 
with respect to the application of our formulas to the calculation of 
the air forces for vibrating lifting surfaces are ommitted since these 
calculations would be essentially identical with those by Schwarz. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 


xs 


L. Schwarz, elsewhere. 
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